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Abstract. 

The dark matter distribution on small scales may depend on the the properties of the first generation of dark matter halos to 
form, which is in turn determined by the microphysics of the dark matter particles. We overview the microphysics of WIMPs 
and our calculations of the collisional damping and free streaming scales. We then plot the resulting density perturbation power 
spectrum and red-shift at which typical halos form, taking into account the effect of uncertainties in the WIMP properties (mass 
and interaction channel) and the primordial power spectrum. Finally we review recent developments regarding the properties 
and fate of the first WIMPy halos. 
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WHY? 

Diverse cosmological observations indicate that the Uni- 
verse contains a significant amount of non-baryonic cold 
dark matter (CDM). Weakly interacting massive particles 
(WIMPs) are an attractive CDM candidate, as a stable 
relic from the electroweak scale genetically has an inter- 
esting present day density, fi c d m ~ ). 

According to standard lore the distribution of CDM is 
independent of its nature (WIMPs, axions, or something 
more exotic). This is true on large (super-kpc) scales, 
but may not be the case on small (sub-pc) scales. Fur- 
thermore the signals expected in WIMP detection exper- 
iments depend on the dark matter distribution. WIMP di- 
rect detection experiments (which aim to detect WIMPs 
in the lab via their elastic scattering off detector nuclei) 
probe the dark matter distribution on sub-milli-pc scales. 
WIMP indirect detection involves searching for the prod- 
ucts (y-rays, antiprotons and neutrinos) of WIMP anni- 
hilation and, since this is a two-body process, the event 
rates are proportional to the density squared and are 
hence enhanced by substructure. Reliable predictions for 
the expected signals therefore require an understanding 
of the distribution of dark matter on small (sub-galactic) 
scales. 

In CDM cosmologies structure forms hierarchically; 
small halos form first with larger halos forming via merg- 
ers and accretion. The density perturbations on very 
small scales, and hence the properties of the first gener- 
ation of structures to form, depend on the micro-physics 
of the CDM and the present day density distribution may 
retain traces of these first structures. 



WIMP MICRO-PHYSICS 

In the early Universe WIMPs are kept in thermal equi- 
librium with the radiation component via chemical inter- 
actions. As the Universe expands the WIMP density de- 
creases, these interactions cease and the WIMPs chem- 
ically decouple or 'freeze-out'. After this the WIMP 
comoving number density remains constant, however 
WIMPs continue to interact kinetically, via elastic scat- 
tering, with the radiation component and are kept in 
local thermal equilibrium. Eventually the WIMP den- 
sity becomes so low that elastic scattering ceases and 
the WIMPs kinetically decouple [l|, |2|, |3J. The aver- 
age momentum exchanged per collision is small how- 
ever so the relaxation timescale, which characterises the 
time at which WIMPs kinetically decouple, is signifi- 
cantly larger than the elastic scattering timescale ulyD. 
The kinetic decoupling temperature was first estimated in 
Ref . yl while Ref . 1 2 ] calculated the last scattering tem- 
perature for supersymmetric WIMPs. While chemical 
decoupling happens at a temperature r c( j ~ GeV), 
kinetic decoupling is delayed by the large entropy of the 
hot Universe and takes place at a temperature T^d ~ <^(10 
MeV). This is generic for any WIMP that is of cosmolog- 
ical relevance today. 

Elastic scattering processes around lead to viscous 
coupling of the WIMP and radiation fluids, resulting in 
collisional damping of WIMP perturbations due to bulk 
and shear viscosity |3J]. The collisional damping scale 
was first derived in Ref. 1 3] by two of the present authors, 
and we presented a more intuitive, but less rigorous, 
derivation using the linearised Navier-Stokes equation in 
Ref. 0| . The net result is exponential damping with a 



characteristic wave-number 



POWER SPECTRUM AND Z NL 
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where m is the WIMP mass. 

One caveat is that our treatment is restricted to sub- 
horizon scales. Near the horizon there are additional ef- 
fects, which have recently been studied numerically by 
Loeb and Zaldarriaga |5]. Their results confirm that the 
essential features of the subhorizon damping are cap- 
tured by our calculation, but the additional terms play an 
important role in determining the precise position of the 
maximum of the CDM power spectrum. We agree with 
their conclusion that an accurate calculation (better than 
10% accuracy) must be based on a numerical computa- 
tion including all of the terms. 

After kinetic decoupling the WIMPs free-stream and 
their distribution function is governed by the collision- 
less Boltzmann equation. We study the subsequent addi- 
tional damping by considering small perturbations away 
from the local thermal equilibrium solution, taking into 
account the perturbations present at kinetic decoupling. 
The characteristic free-streaming scale becomes constant 
after matter-radiation equality and is given by 

_ 1.70x 10 6 (m/100GeV) 1 / 2 (7k d /30MeV) 1 / 2 
fs ~~ Mpc l+ln(r kd /30MeV)/19.2 

(2) 

The free streaming leads to exponential damping of the 
WIMP density contrast again but with an additional poly- 
nomial pre-factor [6, 4]. Since kf s < k^, the cut-off in the 
power spectrum is determined by the free streaming scale 

kf s . 

We focus on four benchmark models which span 
the range of most plausible WIMP properties and have 
present day densities in accordance with the WMAP 
measurement of the cold dark matter density. The de- 
tails of these benchmark models, including the values 
of kd and kf s , are tabulated in Table 1. I parameterises 
the temperature dependence of the thermally averaged 
elastic scattering cross section: ((7 e i) = <7q(T /m) l+l . In 
the Standard Model, elastic scattering between a light 
fermion and a heavy fermion is mediated by Z° exchange 
and / = 0. Other channels may occur however. In super- 
symmetric extensions of the Standard Model, where the 
lightest neutralino is a WIMP candidate, sfermion ex- 
change occurs (and if the neutralino is a gaugino, Z° ex- 
change is suppressed), in which case 1 = 1. Models B and 
C are very close to the bino-like neutralino case consid- 
ered in Ref. |6]. Models A and D show that there is more 
spread in the predicted damping scale if the assumption 
of a bino-like WIMP is dropped. 



The (processed) power spectrum depends on two more 
ingredients: the gravitational growth of perturbations and 
the primordial power spectrum. We calculate the trans- 
fer function, which encodes the gravitational growth, 
by solving the equations which govern the evolution of 
perturbations in matter and radiation fluids during two 
overlapping epochs, radiation domination and close to 
matter-radiation equality, and then matching these solu- 
tions together (see Ref. |4j for further details). We ne- 
glect perturbations in the baryon fluid. This approxima- 
tion is valid on small scales k > k\, ~ 10 3 Mpc _1 for 
z > Zb ~ 150. We have verified that our analytic expres- 
sions are accurate at the 10% level by comparing them 
with the output of the COSMICS package. This is suffi- 
cient given current uncertainties in the cosmological pa- 
rameters and the primordial power spectrum. 

The resulting linear power spectrum, defined as 
^A(k,z) = k 3 /(27t 2 )(\A(k,z)\ 2 ), at z = 300 (before 
the onset of non-linear structure formation) is plotted 
in the left panel in Fig.l for our benchmark WIMP 
models and a WMAP normalised scale-invariant primor- 
dial power spectrum. The characteristic free-streaming 
wave-number is smallest for model A (Dirac WIMP 
with 1=0) and increases with increasing mass for the 
Majorana WIMPs (models B-D). This is clearly reflected 
in the position of the cut-off in the power spectrum in 
Fig. 1. 

The cut-off in the WIMP power spectrum sets the scale 
of the first halos to form. We estimate the redshift at 
which typical (l-ff) fluctuations on comoving scale R go 
nonlinear via o(R,z n \) = 1, where <j(R,z) is the mass 
variance defined by 

r°° (\k 
a 2 (R,z) = / W 2 (kR)^ A (k,z)-, (3) 
Jo k 

where W(kR) is the Fourier transform of the window 
function, which we take to be a top hat, divided by its 
volume. 

The right panel of Fig. 1 shows Zn\ as a function of 
scale R. The cut-off in the processed power spectrum at 
k — 10 6 Mpc _1 leads to a plateau with z nl = zjjf* at R < 
R min = ) p C . We emphasise that z™ ax is the redshift at 
which hierarchical structure formation starts at a typical 
place in the universe, the very first WIMPy halos will 
form significantly earlier from rare large fluctuations. 
The order of magnitude variation in kf s leads to a similar 
variation in /? m ; n and also (because of the dependence of 
the amplitude of the peak of the power spectrum on the 
cut-off scale) a range of values z™ ax « 50 to 65. 

We now turn our attention to the effect of the uncer- 
tainty in the primordial power spectrum. On the scales 
probed by the CMB [^(0.01 - O.^Mpc" 1 ] the primor- 
dial power spectrum is close to scale invariant. The free- 



TABLE 1. Benchmark WIMP models. 



Ref. 


I 


m (GeV) 


r cd (GeV) 


Tied (Mev) 


Mpc" 1 ) 


Mpc l ) 


A 





100 


3.6 


1.6 


14 


0.42 


B 


1 


50 


1.9 


21 


39 


0.94 


C 


1 


100 


3.7 


25 


61 


1.5 


D 


1 


500 


17 


37 


180 


4.0 
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FIGURE 1. The processed power spectrum at z = 300, 3?i±{k,z = 300), (left) and the red-shift at which 1 — rj fluctuations go 
non-linear, z n h (right) for a WMAP normalised scale-invariant primordial power spectrum and from left to right/bottom to top a 
Dirac WIMP with m = 100 GeV and Majorana WIMPs with m = 50, 100, 500 GeV. Without the effects of collisional damping and 
free streaming these quantities would be given by the dotted line. The vertical dashed line denotes k\,, the wavenumber below which 
baryons follow CDM. Our approximations are optimised for k > k\,. 
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FIGURE 2. As Fig. 1 for WIMP benchmark model C (a bino like WIMP with m = 100 GeV) and from top to bottom false 
vacuum dominated hybrid inflation (n = 1 .036, a = 0), a scale invariant primordial power spectrum (n — 1 , a = 0), power law 
inflation (n = 0.964, a = 0) and m 2 (f> 2 chaotic inflation (n = 0.964, a = -0.0006). 



streaming scale kf s ~ 10 6 Mpc~' is seven orders of mag- 
nitude smaller and therefore even a very small scale 
dependence of the power spectrum could significantly 
change the amplitude of the power spectrum at the cut- 
off scale, and hence the red-shift at which the first WIMP 
halos form. The most commonly used parameterisation 
of the primordial power spectrum is 

/ j, \ »(*o)-l+!a(*oMV*o) 

= nh) [j-j , (4) 

where a(k) = dn/dk. The spectral index n and its run- 
ning a depend on the inflationary potential. 

To assess the effects of possible scale dependence of 
the primordial power spectrum we consider three bench- 
mark inflation models which span the range of possi- 
ble power spectra for simple inflation models: V — m 2 (j) 2 
chaotic inflation, with N = 55 e-foldings of inflation so 
that n—\ = —0.036 and a = — 6 x 10~ 4 , power law in- 
flation (a « t p ), with p = 55.4 so that re — 1 = -0.036, 
to match the chaotic inflation model, but in this case 
a = and false vacuum dominated hybrid inflation with 
the largest spectral index allowed by the WMAP data, 
re — 1 = 0.036 and a = 0. For further details see Sec. 7 of 
Ref. H. 

The amplitude of the primordial power spectra on the 
free-streaming scale ~ 10 6 Mpc _1 varies by a factor of ~ 
2.5 (equivalently the amplitude of the fluctuations varies 
by ~ v2-5 ~ 1.6). The resulting processed power spectra 
and the red-shift at which scales go non-linear are plotted 
in Fig.2 for a Majorana WIMP with m = lOOGeV. The 
variation in the amplitude of the primordial power spec- 
tra at the cut-off scale translates directly into a variation 
in the peak amplitude of the processed power spectra and 
hence z"{ ax . 



FATE OF THE FIRST HALOS 

We estimate the size and mass of the first generation of 
typical subhalos that form at z™ ax using the spherical col- 
lapse model The mean CDM mass within a sphere 
of comoving radius R is M(R) = 1.6 x lO~ 7 M (fi/pc) 3 , 
with over-densities having twice this mass. The present 
day radius, after turn-around and violent relaxation, 
would be of order tens of milli-pc, comparable to the size 
of the solar system, and smaller. 

Numerical simulations are required to study the de- 
tailed properties, and fate, of the first generation of 
WIMPy halos to form. Diemand et al. ||7|] have recently 
published the results of the first such simulations to be 
carried out. They use a multi-scale technique to study 
the non-linear collapse of a small, (60 comoving pc) 3 , 
average density region starting at z — 350. The first non- 
linear structures form at z ~ 60 and have M ~ 1O~ 6 M , 



in good agreement with our analytic calculations. The 
simulations have to be stopped at z ~ 25 when the high 
resolution region begins to merge with the surrounding 
lower resolution regions. 

These mini-halos are subsequently subject to a num- 
ber of dynamical processes which could destroy them. 
Matter will be tidally stripped from the outer regions if 
the gravitational force of the parent halo exceeds that of 
the sub-halo. Diemand et al. estimate that the mini halos 
should not be tidally destroyed outside the inner few par- 
sees of the Milky Way |7], in contrast to earlier analytic 
calculations by Berezinsky et al. |8]. The mini-halos can 
also be destroyed via interactions with compact objects. 
Zhao et al. argue that encounters with stars destroy all of 
the mini halos in the solar neighbourhood |9], Moore et 
al. however disagree [10]. Further numerical simulations 
of the dynamical evolution and interactions of mini-halos 
are probably required to resolve these issues. 



SUMMARY 

WIMP direct and indirect detection both probe the dark 
matter distribution on small scales, which depends on 
the properties, and hence fate, of the first generation 
of WIMP halos to form. Collisional damping and free- 
streaming erase density perturbations on small scales 
(k > 10 6 Mpc~') and set the scale of the first generation 
of WIMP halos to form (Af ~ 1O" 6 M ). The fate of these 
halos, specifically whether or not a significant fraction 
survive to the present day within the Milky Way, remains 
an open question. 
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